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Abstract 

We give a combinatorial description of the stationary measure for 
a totally asymmetric exclusion process (TASEP) with second class 
particles, on either Z or on the cycle Zat. The measure is the image 
by a simple operation of the uniform measure on some larger finite 
state space. This reveals a combinatorial structure at work behind 
several results on the TASEP with second class particles. 

1 Introduction 

In the normal totally asymmetric simple exclusion process (TASEP), a num- 
ber of particles occupy some distinct vertices of a graph, which in our case 
will be taken to be either Z or the cycle Zjv- Each particle moves to the posi- 
tion to its right at rate one whenever that position is unoccupied. This gives 
a Markov chain with continuous time on the set of configurations — sets of 
locations for the particles. It is well known (and easy to prove) that if there 
are a particles on Zjy, then the stationary distribution is uniform among all 
( N } possible states. Since in the cyclic case the total number of particles is 
invariant, any other stationary measure is some linear combination of these 
uniform measures. Since the set of stationary measures is convex, from now 
on we focus only on the extremal stationary measures. 
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In the infinite case of Z there are two types of extremal stationary mea- 
sures. There are the so called trivial stationary measures where the particles 
occupy all positions to the right of some point and no other positions. The 
non-trivial stationary measures are those where each vertex is occupied with 
some probability p £ [0, 1] independently of all other vertices. Note that 
these measures are limits as iV — > oo of the unique stationary measure on 
Zjy with \pN] particles. For additional background on the TASEP and many 
references see [HJ [7] . 

We are interested in the process with second class particles defined as 
follows (see e.g. |2j). The particles are classified as first class and second 
class particles. Sites are in one of three states: empty, or occupied by a single 
particle of one class or the other. As with the simple process, each particle 
jumps to the position to its right at rate 1 when that position is empty. 
Additionally, whenever a first class particle has a second class particle to its 
right, the two swap places at rate 1 (thus second class particles may move in 
both directions). 

It is interesting at this point to note that the TASEP on a cycle is equiv- 
alent to the following shuffling method of cards. The cards are arranged in a 
cycle. Each consecutive pair is chosen at rate 1, and the pair is sorted with 
the larger card to the right. Assume there are a high cards, b medium cards 
and c low cards, and that cards within the same class are not distinguished 
from one another. The dynamics are exactly equivalent to those of the ex- 
clusion process with high cards corresponding to first class particles, medium 
cards to second class particles and low cards to empty positions (see [T] for 
an application). It is of course interesting to study the case where there are 
more than three types of cards (or particles). These cases seem significantly 
harder to analyze, though experimental studies of small well as the 

extremal case of iV different particles do show interesting phenomena. 

We give one final equivalent interpretation of the process, that may be 
related to the fact that adding further classes of particles breaks known 
techniques for working with the model. Non-empty sites of the graph are 
occupied by either a particle or an anti-particle. Each edge is chosen at rate 
1. If a particle can move right across the edge to an empty spot it does. If an 
anti-particle can move left across the edge to an empty space it does. Finally, 
if there is a particle on the left and an anti-particle on the right, they both 
move, exchanging their positions, and otherwise nothing happens. Thus there 
are particles moving right and anti-particles moving left, but the movements 
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have rate 1 at each edge rather than each particle. This is equivalent to 
the previous form, with anti-particles representing empty spaces, and empty 
spaces representing the second class particles. When writing out states of the 
process, we will use l's for particles, O's for anti-particles and *'s for empty 
spaces. We will mostly refer to this interpretation of the process from now 
on, and use the term TASEP to refer to this process as well. 

If the initial state is such that only two of the three types of positions 
appear (i.e. if there are no particles of one type or the other, or if all sites are 
occupied), then the process is simply the well understood exclusion process. 
However, the relation to the 1-type case is deeper than that. If an observer 
sees only the particles, and ignores the anti-particles, treating them as empty 
spaces, then he observes a regular exclusion process. Thus in the stationary 
distribution with a particles and b anti-particles, the marginal of the positions 
of particles is uniform over all ( N ) sets. Similarly, the anti-particles form an 
exclusion process on their own (up to a reversal of the directions) and are 
therefore uniformly distributed over (^) possible sets. Of course the two sets 
of locations are disjoint and hence dependent. 

A similar statement holds in the infinite space case. In a stationary mea- 
sure, if one type of particles are in a trivial state, then so must be the other, 
and we get the trivial stationary measures where a single segment of empty 
spaces has particles on its right and anti-particles on its left. Otherwise, the 
set of positions of particles have an i.i.d. product measure with some density, 
as do the anti-particles. Of course, the two marginals are not independent, 
nor is their joint distribution a product measure. 

The main result of this paper is a combinatorial description of the sta- 
tionary measures for the TASEP with three types of positions (empty and 
two types of particles). These results are then used to shed light on some 
known properties of the stationary measures. Duchi and Schaeffer [5] have 
found a similar relation, and use it to derive other results on the TASEP 
(primarily on a finite interval). 

2 Definitions and Results 

Definition 2.1. Two sets of positions S,T, of positions on Z or Zjv, (not 
necessarily disjoint), are said to collapse to a state x of the TASEP if x is 
the result of the following (collapsing) procedure: Anti-particles are placed 
at the locations specified by T. Next, the locations in S are checked (in an 
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arbitrary order). If a location is empty, a particle is placed there. Otherwise 
a particle is placed in the nearest empty position to the left of the specified 
location. 

In the case of the cycle Zjv, we will only use the process for sets S, T with 
1 51 + |^| — n i so that there is always an open position for every particle in 
S. In the case of Z, if there are no open positions to the left of some element 
in S, we disregard that element. 

The collapsing procedure defines a function from the product space (pairs 
of sets) to the space of states of the TASEP. The order at which the positions 
of S are used has no effect on the final resulting state x: The anti-particles in 
x are located exactly in the positions in T. A position a contains a particle in 
the resulting state x if and only if a ^ T and there is some interval / = [a, b] 
such that | J D S\ + 1 I fl T| > |/|. These statements may be used as an 
equivalent condition of collapsing, avoiding the possibly infinite algorithmic 
definition in the case of Z. 

Theorem 2.2. The stationary measure for the exclusion process on Z^r with 
a particles and b anti-particles is the image by collapsing of the uniform 
measure on pairs of sub-sets S, T of the cycle of sizes a, b respectively. 

As a corollary, it follows that the least likely states for given population 
sizes are 1 . . . 1* . . . *0 . . . and its cyclic shifts, i.e. positions where sites of 
each type form a single interval, with positive particles to the left of the 
empty sites (and anti-particles to the right). These states have probability 
1/ (^) (T) each, while every other state has a probability that is some integer 
multiple of this probability. Section El contains some conjectures generalizing 
these facts to the case of more particle types. 

As noted above, On Z the TASEP with second class particles has trivial 
stationary measures, which correspond to the states with particles to right 
of some point, and anti-particles to the left of some other point with empty 
spaces in between. In the notation of first and second class particles these 
are states with second class particles in some interval, first class particles to 
the right and empty spaces to the left. 

Theorem 2.3. The non-trivial extremal invariant measures for the exclusion 
process on Z are the image by collapsing of sets S, T where each n G S with 
probability p and n G T with probability q, all independently. 
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Note that if p + q > 1, then a.s. all positions of Z will eventually have 
some particle on them, and then the resulting measure is one where each 
site has an particle with probability 1 — q, and an anti-particle otherwise, 
independently of all other sites. As noted above, in the lack of empty spaces 
these are the (non-trivial extremal) invariant measures. If p+q < 1, then the 
resulting measure has particles with density p, anti-particles with density q 
and empty sites with density 1 — p — q. 

Theorems 12.21 and 12.31 can also be used to derive properties of the col- 
lapsing procedure. Since the distribution of the positions of particles in the 
stationary distribution is known: 

Corollary 2.4. Let S' is the set of locations of particles after collapsing 
independent S, T , which are have either the uniform distribution over sets of 
a given size (in or i.i.d. product distributions (on Then S' also has 
a uniform or an i.i.d. product distribution. 

Furthermore, the TASEP itself has a natural symmetry of reversing the 
charge of the particles as well as the direction of the graph. It follows that a 
similar "dual" collapsing procedure where the particles are fixed and the anti- 
particles are moved forward to empty spaces also yields the same stationary 
measure. 

In the next section we prove a combinatorial lemma that is closely related 
to stationarity of the collapsed uniform measures. Sections 0] and contains 
the proof of Theorems 12.21 and 1231 respectively. In Section |5] we also use the 
collapsing description of the stationary measures to shed new light on some 
of the results of 0]. Finally, Section |U] contains some open problems and con- 
jectures regarding more general multi-type asymmetric exclusion processes. 

3 Binary trees and dominated sequences 

This section contains the combinatorial foundation for proving the station- 
arity of the collapsed uniform measure. The key result here is a bijection 
between binary trees and pairs of binary sequences, (i.e. sequences made 
up of O's and l's). The sequences are related to the TASEP since a binary 
sequence describes a segment in a state of the TASEP with no empty sites 
(which is one of the reasons we use the {0,*,l} notation). 

Definition 3.1. Consider two finite binary sequences A, B of the same length 
n. We say that A dominates B and write A !>= B if it is possible to get from 
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A to B by moving l's to the right. The weight of a binary sequence A is 
defined as the number of binary sequences dominated by it: 

W(A) = \{B : A )p B}\. 

In particular, it is necessary for A )p B that both sequences have the 
same number of ones. Let the number of l's in the first % digits of A (resp. 
B) be denoted by a(i) (resp. b(i)). A condition equivalent to A )p B is that 
a(i) > b(i) for all i, and a(n) = b(n). Thus for example, W(1010) = 5. If 
x = 1 ... 10. . .0 has k ones followed by / zeroes, then W(x) = ( k ^ 1 )- If the 
l's and 0's were exchanged W(x) would have been 1. 

We use the following version of binary trees: 

Definition 3.2. A binary tree is a rooted tree where each vertex including 
the root may have a left child, marked as left, and may have a right child, 
marked as right. A vertex may have either child, both, or neither. 

Thus a tree is either the empty tree with only a root vertex, or it has a 
left sub-tree, a right sub-tree or both. Note that having an empty sub-tree 
on some side is different from not having a sub-tree on that side. 

Next, we define recursively a function / mapping binary trees to binary 
sequences, as follows. The empty tree (with no edges) encodes the empty 
sequence. Otherwise, 

{f(L)0 if T has only a left sub-tree L, 

lf(R) if T has only a right sub-tree R, 

f(L)01f(R) if T has sub-trees L, R, 

where e.g. f(L)01f(R) means a concatenation of f(L), 01 and f(R). In this 
way any binary sequence may be encoded by a binary tree, though generally 
the encoding is not unique. The length of the sequence is the number of 
edges of the tree, and the number of l's is the total number of right children. 
The following combinatorial Lemma and immediate Corollary relate binary 
trees and dominated sequences: 

Lemma 3.3. There exists a bijection between binary trees T and pairs A, B 
of binary sequences such that A !j= B and A = f(T). 

Corollary 3.4. The number of binary trees that encode a given binary se- 
quence A is W(A). 
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Proof of Lemma Lemma \3.<A First, we define in a similar fashion a second 
function g, that maps binary trees to binary sequences. The empty tree is 
mapped to the empty sequence, and otherwise, 



We now show that mapping a tree T to a pair of sequences A = f(T), 
B = g(T) gives a bijection satisfying the Lemma's requirements. 

The following facts are clear: for sequences A,B,C,D such that A ^ B 
and C ^= D we have AO )p OB and AC )p BD. By induction, it follows 
that f(T) ^ g(T): If the root has only a right child, then f(T) = lf(R) )p 
lg(R) = g{T). If the root has only a left child, then f{T) = f(L)0 ^ Og(L) = 
g(T). If the root has two offspring, then concatenation of the previous two 
cases yields f(T) )p g(T). 

To see that this mapping is a bijection, we show how to (recursively) 
recover from sequences A,B a tree that is mapped to them, in such a way 
that at each step there is a unique possibility, so that the tree is unique. 

Consider a binary tree T, with two sub-trees L, R and sequences A = 
f(T) = f(L)01f(R) and B = g(T) = Og(L)lg(R). Recall that a(i) and 
b{i) count l's in preambles of A and B. Since f{L) )p g(L), we have that 
o(0 — b{i + 1) for i < \L\. However, a(\L\ + 1) = b(\L\ + 1) is the number 
of right edges in L, and the next bit in B is a one. Thus the inequality 
o(0 — b(i + 1) fails for for the first time for i = \L\ + 1. The following 
algorithm emerges: Given A, B find the first i for which a(i) < b(i + 1), 
and set \L\ = i — 1. This generally identifies a unique representation of 
the sequences as A = X01Y and B = OX'lY' with X ^ X' and Y ^ Y' . 
To reconstruct the unique tree T mapped to (A, B), proceed recursively to 
identify L from X, X' and and R from Y, Y' . 

It remains to see that the cases where the above procedure fails to locate 
a representation as above correspond exactly to cases where T has only a 
left or only a right sub-tree. One possibility is that a(0) = < 1 = b(l), 
and then the above would give \L\ = —1. In this case A = 1Y and B = 1Y' 
with Y )p Y', so the unique tree mapped to A, B has only a right sub-tree 
R, where R is mapped to Y, Y'. 

The other extreme is the case that a(i) > b(i + 1) for all i < n, suggesting 
\R\ = —1. In this case we have A = XO and B = OX' with X )p X', and we 




Og(L) if T has only a left sub-tree L, 

lg(R) if T has only a right sub-tree R, 

Og(L)lg(R) if T has sub-trees L, R, 
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find a tree with only a left sub-tree. Thus in all cases the algorithm proceeds 
recursively to find a unique tree that is mapped to A, B. □ 

Lemma 3.5. For a binary sequence A, 

W{A) = W(X)l A=X0 + W{Y)1 A=1Y + W(X)W(Y). 

X01Y=A 

Thus if A ends with a and equals XO, the RHS gets a contribution 
of W(X) from the first term. Similarly, if A begins with a 1 there is a 
contribution from the second term. The sum in the RHS has a term in the 
sum for each representation of A as X01Y. For example, if A = 1011010 
then we get 

w(A) = w(ionoi) + w(onoio) + w(i)w(ioio) + w(ion)W(p) 

= 7 + 9 + 1-5 + 2-1 =23. 

Proof. By Corollarv l3.4l it suffices to show that the RHS equals to the number 
of binary trees that encode A. This is done by induction. If A = XO ends 
with a 0, then there are W(X) trees encoding A with only a left sub-tree 
(which can be any tree encoding X). If A = 1Y then similarly there are VK(Y") 
trees encoding A with only a right sub-tree. Finally, for any occurrence of 
01 in A where A = X01Y there are W(X) possible left sub-trees and W(Y") 
possible right sub-trees, giving a term in the sum. □ 

4 Proof of Theorem 12.2 

Consider a state x of the exclusion process on the cycle. How many pairs 
of sets S, T collapse to x? Since the collapsing procedure begins by placing 
the anti-particles at T, the unique T is given by the set of positions of anti- 
particles in x. There may be a number of different sets S that (together 
with T) collapse to the state x. In order for the collapsing process to reach 
x it is necessary that S contains none of the empty positions of x (positions 
marked with *'s). The empty positions in x break up the cycle into a number 
of segments each containing a sequence of particles (l's) and anti-particles 
(0's). Denote the binary segments of x by A\, . . . , A\. 

During the collapsing procedure, if an element p G S results in a particle 
being placed in some position q to the left of p, there can be no empty position 
in the interval [q, p] , since otherwise the particle would have been placed there 
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instead. Thus the elements of S in each such binary segment must collapse 
into the positions marked for particles in that segment. It follows that for 
each binary segment Ai, the sequence having l's at the elements of S in that 
segment is dominated by Ai, and so there are W(-Ai) possibilities for the 
intersection of S with that segment. The total number of possibilities for S 
is therefore Y\ W(Ai), and the collapsed uniform measure of the state x is 

P(x) — /jy\ 7JK 



For example the cyclic state * 10** 10 100*0 101 may be reached from 

W{l0)W{(f))W{l0100)W{010i) = 2 • 1 ■ 9 ■ 2 = 36 

sets S and so its probability is 36/ Q 5 ) (g 5 ). 

To show that the collapsed uniform measure is stationary, place at each 
state x a mass m(x) = f] W(Ai) - - a multiple of the collapsed uniform 
measure. Let the mass flow according to the transition kernel of the process 
(so if the process passes from x to y at rate r, mass flows from x to y at a 
rate of r ■ m(x). It suffices to show that the derivative of the mass at any 
state x is 0. 

Denote x — > e y if an action (sorting particles) along an edge e leads from 
state x to state y. Since mass flowing from x to itself makes no difference, 
we only use this notation for x ^ y. Since the edge is determined uniquely 
by x and y, the edge subscript will usually be omitted. Each edge is used at 
rate 1, so the mass derivative is given by 

jm(x) = m (y) - = m ^ 2 ( ~ Yl 1 ) • w 

We now associate each term in each of the sums of (0) with one of the 
binary sequences Ai appearing in x in such a way that the terms associ- 
ated with each sequence will cancel out, proving the Theorem. Since edges 
connecting two empty positions lead from x to itself and have been disre- 
garded, each edge intersects exactly one of the binary segments of x, and the 
corresponding term is associated with that sequence. 

The mass in the first sum — flowing into x — corresponds to edges with 
end-points marked '01', '*1', or '0*' in x. For each such edge we need to 
calculate the mass at the state resulting from 'unsorting' the edge. When 
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such an edge is unsorted, the resulting state y is very similar to x. Indeed, 
{Bi} are the binary sequences appearing in y, then all but at one or two of 
them are equal to those in x. 

• Consider first the case where y — > e x and the endpoints of e are marked 
'0*' in x and '*0' in y. In this case y has the same binary sequences as x 
except for two: A { = B { and B i+1 = 0A i+1 . Since W(0A) = W(A), it 
follows that = 7tftt4, where Bj is A- with a terminating removed. 

m(x) W(A i ) ) % % o 

• Similarly, if y —> e x and the endpoints of e are marked '*1' in a; and '1*' 
in y, In this case A, = 1B { and B^ = A^l. Since W(Al) = W(A), 
we find that = ^[ B A , where Bi is A; with an initial 1 removed. 

m(x) W (Ai) ' 1 1 

• Finally, if e is marked with '01' in x, then spliting around that edge 
we have Ai = X01Y and Bi = X10Y, with all other sequences beeing 
equal. Again, = fgj. 

Let us extend the — > notation to binary sequences, so that A — > B if it 
is possible to pass from A to B by either removing an initial 0, removing a 
terminating 1, or replacing a 10 by 01 somewhere in A. Consider the terms 
in the RHS of (0) that are associated with the binary sequence Ai of x. After 
substituting the above for ^r4, these terms come to 

° m(x) ' 

W(B) 



mix 



E 




and so it suffices to prove for an arbitrary sequence A that 

W{B) = W(A). (2) 

B^A A^C 

Given a binary sequence A, associated terms in the RHS correspond to 
each occurrence of '10' in the sequence, as well as to an initial if there is 
one and to terminating 1 if the sequence ends with a 1. The number of such 
terms is always one more than the number of times '01' appears in A. 

The terms in the LHS, are W(B) where B results from A either by replac- 
ing '01' by '10' at some place, or by removing an initial 1 or a terminating 
if A has them. Consider a pair of sequences B —>■ A where A = X01Y and 
B = X10Y. Since B )p A, B also dominates any sequence that A dominates. 
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The difference W(B) — W(A) corresponds to sequences that B dominates and 
A does not. Such a sequence must be of the form X'lOY' , where X ^= X' 
and Y )? Y', and so W(B) = W(A) + W(X)W(Y). 

Substituting this for of W(B) in (J2J) results in cancellation of all but one 
term on the RHS, and the resulting needed identity is exactly that given by 
Lemma 13.51 

5 The infinite setting 

Proof of Theorem \2.?A It is known (see jS] ) that there is a unique non-trivial 
stationary measure with marginals p, q for particles and anti-particles, and 
that those are all the non-trivial extremal stationary measures. It only needs 
to be shown that the collapsed i.i.d. measures are stationary and have the 
correct marginals. 

The collapsed i.i.d. measure is the limit as N — > oo of the collapsed 
uniform measure on a cycle of length N with \pN] particles and [qN] anti- 
particles. Since in the finite case particles and anti-particles are uniformly 
distributed over all subsets of the appropriate size, in the limit they have 
densities p and q respectively. Since correlations in this measure decay ex- 
ponentially in the distance, it follows that the limit is stationary. □ 

The collapsing procedure sheds new light some known results. First, the 
fact that second class particles induce a factoring of the stationary distribu- 
tion (see 0): given that is an empty position, the state in Z + and in Z~ 
are independent. With the collapsing procedure, the state in Z + depends 
only on the positive elements of S, T. Conditioned on the event that re- 
mains empty, no particle crosses from Z + to Z - , and so the state on Z~ is 
determined by the negative elements of S, T and is therefore independent of 
the state on Z + . 

Next, consider the relation demonstrated in [HI Ej to a certain biased 
random walk on Z. For marginals p for particles and q for anti-particles, the 
corresponding random walk has steps of —1, 0, 1 with step distribution given 
by P(X = 1) = (1 - p)(l - q), and P(X = -1) = pq. When p + q < 1 we 
have EX = 1— p — q > 0. 

The random walk is very naturally coupled with the pair S,T of i.i.d. 
subsets of Z with densities p, q, and is given by Z n = \T C R [—n, 0)| — \S fl 
[— n, 0)| (note that the sets are explored backwards), which clearly has the 
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above step distribution. Lemma 2.5 of [3] says that the distribution of the 
distance between second class particles is the same as the hitting time of 1 by 
the random walk. In our notation this is a statement regarding the distance 
between empty positions. Given S, T such that there is no particle at 0, the 
next hole to the left of is at — n exactly when the corresponding Z hits 1 
at time n. 

The following is standard and easily seen: The stationary measures seen 
from a single second class particles are derived from the stationary measures 
with anti-particles, by the following rule. Condition on having an empty 
position at 0, and place a second class particle there. All particles are first 
class particles, and any empty spot to the right of is also filled with a first 
class particle. All other positions become empty. By selecting the densities 
of particles and anti-particles the asymptotic densities to either side can be 
controlled. 

Extend the definition of the random walk to negative indices by Z_ n = 
— |T c n[0, , n,)| + |5'n[0,n)|. Theorem 5 of |I] states that the stationary measure 
for the process seen from a single second class particle with given asymptotic 
densities is very close to a product of independent i.i.d. measures on Z + 
and Z~ with the given densities. Formally, the stationary measure and the 
product measure can be coupled so that they are exponentially unlikely to 
differ in many positions. 

This again follows from the collapsing procedure, since the final location 
of a particle is exponentially unlikely to be far from the position given by S. 

6 Open problems 

The following conjectures are a generalization of our results to more classes 
of particles, both in qualitative and quantitative terms. 

Conjecture 1. In a cycle with particles of a number of classes, in the sta- 
tionary distribution, the least likely states are those where if the cycle is cut 
at some point the particles are arranged in reversed order of speed. 

Let x be one of the above states (there are N of them). 

Conjecture 2. In the stationary distribution, P(x) = Y[ ( N .) > where n is 
the cycle length and Si is the number of particles of class at least i. 
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Conjecture 3. In the stationary distribution, the probability of any other 
state is an integral multiple of P(x). 

Finally, is there a useful generalization of the collapsing procedure for 
processes with more classes of particles? The obvious generalizations of re- 
peated collapsing do not appear to give the correct stationary distribution. 
The smallest case where they fail is that of a cycle of length 4 with all 
different particles (or 3 particles and an empty space). It turns out that 
//(1324) 7^ /i(1423) under the stationary measure /i. 
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